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ABSTRACT 


A  least-squares  method  for  evaluating  aircraft  stability 
derivatives  and.  transfer  function  constants  from  flight  test  fre¬ 
quency  response  data  is  derived.  The  normalized  least-squares 
equations  for  rotating  vectors  are  derived  by  minimizing  the  sum  of 
the  squares  of  the  residual  vectors  without  regard  to  the  phase 
angles  of  these  residuals. 

This  method  was  applied  to  longitudinal  oscillation  data  ob¬ 
tained  from  flight  tests  of  a  USAF  bomber-type  B-25J  aircraft.  It 
is  shown  that  a  linear  dependency  in  the  longitudinal  variables  of 
motion  must  be  resolved  before  any  solution  can  be  obtained.  The 
assumption  of  a  linear  relationship  between  the  stability  derivatives 
for  the  forces  or  moments  resulting  from  the  rate  of  change  of  the 
angle  of  attack  and  the  pitching  velocity  eliminated  this  dependency 
and  gave  acceptable  results.  The  lift  caused  by  pitching  velocity 
must  be  included  in  order  to  obtain  a  reasonable  value  for  the  lift 
resulting  from  elevator  deflection. 

Close  agreement  was  obtained  for  the  related  values  of  the 
transfer  function  constants  from  the  data  of  several  independent 
measurements. 

The  application  of  vector  algebra  greatly  simplifies  the  numerical 
computation  required  to  obtain  the  normalized  equations. 
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INTRODUCTION 


Considerable  effort  has  been  expended  in  recent  years  to  improve  flight 
test  methods  of  evaluating  aircraft  stability  derivatives  and  transfer 
functions.  These  improved  methods  are  required  to  provide  more  accurate  data 
for  the  design  of  automatic  control  equipment  and  to  give  the  aircraft  de¬ 
signer  a  better  check  on  wind  tunnel  and  estimated  parameters „  Satisfactory 
aircraft  handling  qualities  are  also  usually  expressed  in  terms  of  some 
stability  parameter,  thereby  necessitating  flight  tests  to  determine  compli¬ 
ance  with  specifications  * 

Most  flight-test  methods  of  evaluating  aircraft  stability  parameters  are 
primarily  steady-state  maneuvers*  A  standard  approach  is  to  maintain  constant 
as  many  variables  of  the  aircraft  motion  as  possible  while  measuring  the 
effect  of  the  remaining  variables*  For  example,  the  slope  of  the  lift  curve 
is  usually  found  by  measuring  the  change  in  the  angle  of  inclination  of  longi¬ 
tudinal  axis  with  changes  in  the  air  speed  in  straight  and  level  flight*  This 
method  is  subject  to  various  possible  errors,  however*  Lift  resulting  from  the 
change  in  elevator  position  and  lift  caused  by  changes  in  the  flow  field  re¬ 
sulting  from  propeller  wwashn  and  airspeed  changes  are  included  in  addition  to 
lift  resulting  from  change  in  angle  of  attack*  Also,  when  the  lift  coefficient 
is  a  function  of  the  Mach  number,  this  method  cannot  be  applied*  The  elevator 
effectiveness  is  usually  evaluated  in  static  methods  by  measuring  the  elevator 
deflection  required  to  offset  a  known  moment  applied  to  the  aircraft*  This 
moment  change  is  obtained  by  shifting  the  center  of  gravity  location*  Again 
some  error  is  involved  in  this  method,  since  some  elevator  deflection  is  re¬ 
quired  to  overcome  the  pitching  moment  resulting  from  changes  in  the  aircraft 
angle  of  attack*  This  change  in  the  angle  of  attack  is  a  by-product  of  the  load 
shift  between  the  wing  and  the  tail  as  the  position  of  the  center  of  gravity  is 
changed*  These  methods  have  the  further  disadvantage  of  requiring  extensive 
flight- test  time  and  expert  piloting  technique* 

The  method  described  in  this  report  adopts  the  approach  of  activating  all 
of  the  variables  of  the  aircraft  motion  simultaneously  while  retaining  the  oper¬ 
ating  conditions  as  constant  as  possible*  Values  of  the  stability  parameters 
are  thai  deduced  in  accordance  with  the  theory  of  probability*  It  is  necessary 
to  use  an  equation  describing  the  relationship  between  the  response  and  the 
stability  parameters*  This  relationship  can  be  expressed  in  the  form  of  equa¬ 
tions  of  motion  or  as  transfer  function  equations* 

The  flight  test  evaluation  of  the  stability  parameters  by  this  method 
requires  that  each  variable  included  in  the  equations  be  measured  directly  or 
be  derived  from  measured  data*  The  evaluation  of  the  transfer  function  con¬ 
stants  in  the  equation  for  the  rate  of  pitch  resulting  from  elevator  deflection 
requires  that  the  first,  second,  and  third  derivatives  of  the  angle  of  pitch  be 
obtained,  as  well  as  the  elevator  angle  and  its  first  derivative.  With  the  best 
instrumentation  techniques  available  to  date,  it  has  been  found  practicable  to 
measure  only  the  first  pitching  derivative  with  an  acceptable  degree  of  accuracy* 
The  remaining  derivatives  mus t ,  therefore,  be  obtained  from  the  measured  rate 
by  analytical,  numerical,  or  graphical  means* 
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Two  possible  techniques  for  exciting  the  variables  of  the  aircraft  motion 
under  constant  flight  conditions  are  the  transient  method  and  the  sinusoidal 
oscillation  method,,  In  the  transient  method  the  aircraft  is  disturbed  from 
a  steady-state  condition  by  a  step  or  pulse  type  motion  of  the  aircraft  con¬ 
trol  surface*  The  oscillation  method  consists  of  obtaining  the  aircraft 
motion  by  a  sinusoidal  variation  of  the  control  surface  about  a  trim  point 
at  various  frequencies* 

Ihe  transient  method  requires  a  minimum  of  flight  time  and  has  the  fur¬ 
ther  advantage  of  being  applicable  to  flight  conditions  other  than  straight 
and  level,  such  as  climbs,  descents,  or  turns*  However,  the  evaluation  of 
the  variables  of  motion  and  $11  required  derivatives  with  sufficient  accuracy 
to  permit  solving  for  the  stability  parameters  does  not  appear  feasible  in 
the  time  plane*  Iq  the  above  example  of  the  transfer  function  constants  for 
the  rate  of  pitch,  6  can  be  measured  directly*  If  the  time  history  of  the 
rate  of  pitch  is  quite  smooth s  the  second  derivative  of  the  pitch  rate  can  be 
obtained  with  fair  accuracy  by  standard  numerical  or  graphical  methods*  How¬ 
ever,  if  the  third  derivative  is  then  obtained  similarly  from  the  second,  the 
inaccuracies  are  no  longer  tolerable*  For  those  elastic  aircraft  having  bend¬ 
ing  modes  that  appear  on  the  record  in  addition  to  basic  airframe  response, 
even  the  first  derivative  that  is  obtained  by  graphical  or  numerical  methods  is 
of  doubtful  value* 

The  oscillation  method  requires  extensive  flight-test  time  and  is  essen¬ 
tially  restricted  to  straight  and  level- flight  conditions*  However,  the  re¬ 
maining  derivatives  of  the  components  of  motion  can  be  obtained  analytically  ih 
the  frequency  plane  with  ease  and  accuracy  from  the  phase  and  amplitude  of  any 
measured  components*  The  measured  component  can  be  expressed  as  a  rotating 
vector  in  the* complex  plane  by  the  relationship  0  »  |0l  e^o  The  second  deriva¬ 
tive  is  Q  =  =  iu)Q  eia>k„  Consequently,  the  derivative  of  a  rotating  vector 

can  be  obtained  by  advancing  the  phase  angle  by  90  degrees  and  increasing  the 
amplitude  by  the  factor iJ*  Higher  order  derivatives  can  be  evaluated  similarly* 
The  oscillation  method  has  a  further  advantage  in  that  the  recorded  data  can 
be  corrected  for  the  error  introduced  by  the  dynamic  response  characteristics  of 
the  instrumentation  by  applying  phase  and  amplitude  calibrations  at  each  fre¬ 
quency* 

The  aircraft  responses  obtained  by  the  two  methods  above  can  be  related 
analytically  by  Fourier  or  Laplace  Transformation  theory*  This  provides  a 
means  of  taking  advantage  of  the  best  features  of  each  technique*  The  flight 
test  time  can  be  minimized  by  recording  measured  responses  in  the  time  plane 
and  transforming  these  responses  to  the  frequency  plane  to  take  advantage  of 
the  simplified  analyses  procedures  available*  The  transform  theories  assume 
that  the  principle  of  superposition  applies,  that  is,  that  all  components  of 
the  motion  are  additive  and  that  the  system  can  be  represented  by  a  set  of 
linear  differential  equations  with  constant  coefficients*  Therefore,  any  non- 
linearities  existing  in  the  response  characteristics  can,  possibly,  lead  to 
prohibitive  errors  in  the  transformation  from  the  time  plane  to  the  frequency 
plane*  The  validity  of  the  transformation  can  be  checked  by  observing  the 
agreement  in  several  frequency  response  characteristics  obtained  by  analyzing 
transients  having  appreciable  differences  in  frequency  content*  A  numerical 
method  for  accurately  obtaining  the  Fourier  transforms  of  arbitrary  time 
histories  is  given  in  Reference  1* 
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SYMBOLS'  AND  NOTATION 


perturbation  elevator  deflection,  radians 

perturbation  angle  of  pitchy  radians 

pitching  velocity s  d©/dt,  radians /second 

2  2  2 

pitching  acceleration,  d  9/dt  ,  radians /second 

perturbation  angle  of  attack,  radians 

perturbation  flight  path  angle,  radians 

angular  frequency,  radians/second 

air  density,  slugs /feet 

m/©VS,  second 

phase  angle 

downwash  angle  at  the  tail,  radians 

rate  of  change  of  downwash  angle  at  the  tail  with  angle  of  attack,, 
summation  of  scalar  products 
determinant  of  coefficients 
lift  coefficient,  2xag/f  V  S 

rate  of  change  of  lift  coefficient  with  angle  of  attack 

rate  of  change  of  lift  coefficient  with  elevator  deflection 

rate  of  change  of  lift  coefficient  with  pitching  velocity 

rate  of  change  of  lift  coefficient  with  rate  of  change  of  angle  of 

attack 

rate  of  change  of  lift  coefficient  with  angle  of  attack  acceleration 

/  2 

moment  coefficient,  21 y/f  V  SC 

rate  of  change  of  moment  coefficient  with  angle  of  attack 

rate  of  change  of  moment  coefficient  with  elevator  deflection 

rate  of  change  of  moment  coefficient  with  pitching  velocity 
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rate  of  change  of  moment  coefficient  with  rate  of  change  of 
angle  of  attack 


C  rate  of  change  of  moment  coefficient  with  acceleration  of  angle 

^  of  attack 

A  in-phase  componait  (with  subscripts  ec,;  9  0,  n) 

B  out-of-phase  component  (with  subscripts  «cs  4  9  n) 

A,B,C,E  transfer  function  constants  (with  subscripts  0,  1,  2) 

D  differential  operator,  d/dt 

i  imaginary  unit,  T)- 1^  (complex  notation) 

i  unit  vector  along  vertical  axis,  (vector  notation) 

V  ’  true  velocity,  feet/second 

lt  tail  length,  feet 

2 

I  aircraft  moment  of  inertia  about  lateral  axis,  slug-feet 

y 

m  aircraft  mass,  slugs 

2 

a  aircraft  vertical  acceleration,  feet/second 

S  wing  area,  square  feet 

C  mean  aerodynamic  chord,  feet 

n  normal  acceleration,  gps 

g  acceleration  due  to  gravity,  32°2  feet/second'1 


second 
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THEOHT 


Having  obtained  an  accurate  expression  of  the  motion  of  all  variables 
of  interest*  the  problem  consists  of  obtaining  values  for  the  associated 
parameters  of  edch  variable*  namely;  the  stability  derivatives  or  the 
transfer  function  constants.  For  simplicity*  the  theory  will  be  developed 
for  the  stability  derivatives  in  the  lift  equation  assuming  that  lift  will 
result  only  from  the  incremental  changes  in  the  angle-of-attack  and  elevator 
deflection.  Equation  S*  Appendix  I,  will  then  be  written  as 

CLf  *  cl/'  -V  U> 

This  equation  relates  the  lift  forces  to  the  vertical  acceleration  at  any 
instant  of  time*  irrespective  as  to  whether  the  motion  is  any  arbitrary 
transient  or  a  steady-state  sinusoidal  oscillation.  In  the  latter  case*  the 
variables  can  be  considered  for  any  given  frequency  to  be  constant  amplitude 
vectors  rotating  in  the  complex  plane  with  fixed  phase  relationships.  These 
amplitudes  and  phase  angles  can  be  experimentally  determined  for  a  number  of 
different  frequencies  at  a  constant  flight  configuration. 

If  these  data  were  obtained  without  error  of  any  form  and  the  above  equa¬ 
tion  described  the  aerodynamic  and  inertial  forces  completely*  it  would  then 
be  possible  to  solve  for  the  unknown  derivatives  directly.  The  variables 
could  be  resolved  into  components  along  the  real  axis  for  as  many  frequencies 
as  there  are  unknowns,  A  simultaneous  solution  of  these  real  equations  would 
yield  the  accurate  values  of  the  derivatives.  Since  neither  of  these  assump- v 
tions  is  true*  this  approach  leads  to  large  errors  in  practice.  Random 
errors  may  be  present  that  are  larger  than  the  effect  of  some  derivatives  at 
the  particular  time  in  the  cycle  selected  for  taking  real  parts.  The  effect  of 
a  neglected  stability  derivative  will  also  be  large  when  its  associated  vari¬ 
able  has  relatively  large  components  along  the  real  axis  as  compared  to  some 
included  variables.  Standard  curve  fitting  theory  for  scalars  can  be  extended 
to  minimize  these  errors  when  the  variables  are  rotating  vectors. 

The  problan  of  obtaining  values  of  the  stability  derivatives  in  the  sim¬ 
plified  lift  equation  given  above  will  now  be  considered  when  the  phase  angles 
and  amplitudes  have  been  obtained  from  flight  data  for  a  number  of  frequencies. 
Any  values  that  are  found  will  not  satisfy  the  equations  exactly  at  each  fre¬ 
quency  point  at  which  data  are  obtained.  This  is  expressed  algebraically  by 
transposing  all  of  the  terms  in  equation  1  to  the  left-hand  side  and  retaining 
the  residual  resulting  from  the  inequality  on  the  right-hand  side, 

<+■  Ct  +  C_n0  —  v.  *  r  —  1*  3?  ©  o  ,  m 

LJ* I  L  1  1 

where  m  *  the  total  of  frequency  points 
v  -  residual 

Since  all  terms  on  the  left  are  vectors*  the  residual*  v^*  will  also  be  a 
vector  at  each  frequency.  It  is  now  required  to  find  values  of  the  de¬ 
rivatives  that  will  make  the  amplitudes  of  these  residuals  small.  The 


CT 

Tx  1 
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relative  phase  angles  of  these  residuals  that  result  from  minimizing  these 
amplitudes  are  of  no  concern  at  this  point*  Later,  it  will  be  shown  that 
these  phase  angles  can  be  used  as  an  indication  of  the  validity  of  the 
equation  being  fitted* 


It  can  be  shown  (Reference  2)  that  the  most  probable  values  of  the  un¬ 
knowns  will  be  such  as  to  make  the  sum  of  the  squares  of  the  residuals  a 
minimum*  The  square  of  a  residual  vector  is  a  scalar  equal  in  magnitude  to 
the  square  of  the  amplitude  of  the  vector  <> 


o  o  o 


:  V1  +  v2  +  v3  +  m 

'  f(V  cbf  1  E  lf  2* 


(3) 


hi. 


The  condition  that  this  function  be  a  minimum  requires  that  its  partial 
derivatives  with  respect  to  the  unknowns  be  zero* 


•  |ff,  -  2<  Vi  *  V1  +  Vi*!  *  ^ &  *  V2  ^  Cl"2W 

♦  •  •  •  2«yi,  ♦  ctJn  ♦  oLnB)«m  -  0 

or  cj^ac±2  +  CL^SCi^i  *  ■^'cLnio0ci^  1  81  2s  3S  »  ■>  » 

Xk  cf 

Similarly,  when 


00 


(5) 


$£  =  0 
0CT 


Ct  2fxi <f±  +  clF^i2  ~  '^cLni',<c4>  1  ~  2s  5?  - 

GC  cf 


(6) 


m. 


Equations  5  and  6  are  the  normalized  equations  for  the  least-squares  solution 
for  two  unknown  coefficients  of  variable  vectors*  The  indicated  products  of 
the  various  pairs  of  vectors  are  defined  as  the  scalar  or  dot  products 
(Reference  3) »  The  scalar  product  of  two  vectors  is  numerically  equal  to  the 
product  of  their  amplitudes  times  the  cosine  cf  the  angle  between  the  two 
vectors*  This  is  equivalent  to  multiplying  the  vector  of  the  variable  being 
normalized  by  the  respective  components  of  the  other  vectors  along  this 
normalized  vector*  In  this  way  each  individual  equation  is  normalized  for 
each  unknown  parameter  at  that  unique  angle  at  which  the  associated  vector  is 
at  its  peak  or  maximum  value*  It  is  apparent  that  each  parameter  will  have 
the  maximum  effect  on  the  response  when  its  vector  is  at  its  peak  value*  It  is, 
therefore,  logical  that  the  parameter  should  be  so  selected  that  it  yields  a 
minimum  error  at  the  point  where  it  has  maximum  effect*  Basically,  in  the 
equations  given  above,  this  amounts  to  resolving  the  vectors  into  real  parts 
with  the  angle-of-attack  vector  lying  on  the  real  axis  when  normalizing  the 
equation  for  Cr  and  then  resolving  the  vectors  into  real  parts  again  with  the 
elevator-deflection  vector  lying  along  the  real  axis  when  normalizing  the 
equations  for  Cp  * 
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Using  ihe  definition  given  for  the  scalar  products  the  normalized  equations 
can  be  simplifiedo 

|a|  glcos  ^ 

=  Igloos  (S^  -  0^-)  (7) 

=  lad  y|(cos  QOS  +  sin  sin  0^,  ) 

where  !  S  indicates  the  absolute  length  of  the  vector 
is  the  angle  between  these  and  $  vectors 
0pj.  is  the  angle  between  fie  and  the  real  axis 

0^  is  the  angle  between  /  and  the  real  a*is 


Let  Pfer  =  |<od  cos  A^ 

Ba  «  W  sin  y 

Then  oc  •  cf  =  VW  +  §xV 


I  cfl 


cos 


0 


'cf 


!c fl  sin  • 


(g) 


Similarly,  the  remaining  scalar  products  can  be  found  in  terms  of  real  and 
imaginary  components  of  the  complex  variables „  Equations  5  and  6  become 


V 


Cl^  ^±kcf±  +  B‘ 


%  ^i 


^  °i 

h  glI 


C^i 


Vi 


;LV"cTi 


.y 

(9) 

+  Br  B  ) 

(10) 

cf.  n~ 
i  1 

The  amplitudes  of  all  the  vectors  are  usually  given  in  terms  of  a  unit  input,  and 
the  phase  angles  are  all  referred  to  the  input  vector  which,  in  this  case,  is  the 
elevator  deflection*  Therefore,  the  elevator  deflection  vector  will  have  a  unit 


length  lying  along  the  real  axis  in  every  case* 
the  normalized  equations  become 


Then  A, 


CL^(4,  +Ba  )  *  V 


«x. 


T(A  A 
L  «i  n^ 


B  B  ) 
a.  n, 
x  x 


1  and  Br  =  0  and 

(11) 


V 


+  n£T 


*i\s 


i  =  1 


si  2,  3j 


me 


(12) 


The  original  vector  equations  can  also  be  separated  into  real  and  imaginary 
parts  and  these  scalar  equations  normalized  by  the  conventional  least-squares 
method*  When  all  the  real  and  imaginary  equations  are  summed  in  the  normalizing 
process,  the  same  equations  as  given  above  will  result*  This  approach,  however, 
makes  it  difficult  to  visualize  the  basic  validity  and  effectiveness  of  the 
m  etho  d  * 


The  least-squares  method  for  vectors  can  now  be  applied  to  obtain  values 
for  all  of  the  moment  and  lift  derivatives  included  in  the  equations  of 
motion  in  Appendix  X„  It  is  first  necessary  to  resolve  the  linear  dependency 
that  exists  between  the  variables  of  motion,  as  shown  in  the  derivation, 
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before  normalizing  the  equations <>  In  the  simplified  theoretical  evaluation 
of  the  stability  derivatives,  it  is  found  that  the  ec  and  ©  derivatives  are 
directly  related  by  a  constant  equal  to  the  rate  of  change  in  downwash  angle 
at  the  tail  with  change  in  aircraft  angle  of  attack  (Reference  1+) .  This  can 
be  expressed  algebraically  as 


(13) 


\  *  K%  <1W 

d£ 

where  K  =  mr,  the  change  in  downwash  angle  at  the  tail  with  change  in 
aircraft  an&Le  of  attack* 


Substitution  of  equations 
I  will  reduce  the  equations  to 


13  and  llj.  into  equations  7  and  g  of  Appendix 
three  linearly  independent  variables* 


c  c  *  c  -(e  *  Bi) 
■V  m<r  m. 


h&° 


-cLn 


(15) 

(16) 


The  solution  of  these  equations  by  least-squares  will  yield  values  of  the  © 
derivatives  directly*  These  derivatives  will  act  along  the  new  vector, 

(6  -  H*)s  and  thereby  include  the  effects  of  the  «*  derivatives.  Since  theoe 
derivatives  are  relatively  small  the  effects  of  inaccuracies  in  the  assumed 
values  of  K  should  be  minor.  Equations  15  and  16  will  be  written  in  matrix 
notation  for  compactness  before  normalizing. 


-cLn 


•  o 

or  h© 


j  =  m  or  L 


(17) 


The  normalized  equations  for  m  data  points  can  be  obtained  directly  by  scalar 
multiplication  of  each  equation  of  motion  by  each  element  in  the  matrix  of 
variables  in  turn  and  summing  the  results* 


£cc2 
i£  «x 


of 

(©  +  fee) 


(©  +  fex) 

(©  + 


’(©  +  fee) 


£  (6 


(ig) 


r-  1 

•m 

£ -CLn  • « 

*h  ©k 

£-cLn  / 

or 

£h©V 

£-CLn  <•  (Q  +  Kk) 
• 

j  -  L 

g  h©  •  (©  +  fe) 

B 
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Expanding  and  simplifying  the  scalar  products  and  lotting 
«  ls  the  normalized  equations  become  « 

^*iw  + %b©1  i  c.i : 


Ca  a»  + 

<x  © 


A  f  +  iB 


©  ""*6 


+  KiL>  &(Q  +  ife)2 


n  +  W 


—CT  A 
t  n 


5  -h  <A„Ai  *  B„V  +  K<AiAn  *  %®n>  3  -  1 


-fail  Bi 


hK(AcxAe  + 


j  *  m 


Similar  solutions  can  be  made  for  obtaining  values  for  the  transfer 
function  constants  given  in  equations  1 2*  13  and  21*  Appendix  I®  The 
normalized  equations  for  the  constants  in  the  equation  for  Q/^  will  be  de¬ 
rived  first®  Rewriting  equation  12*  Appendix  I,  and  transposing  terms  to 
eliminate  fractions  gives 


AQ©  +  A1D©  +  A2D  § 


=  Bq<£  +  B-^D^ , 


To  obtain  a  unique  solution  for  the  constants*  it  is  first  necessary  to  reduce 
the  number  of  unknowns  by  one.  Since  Ag  is  almost  entirely  made  up  of  known 
quantities*  it  is  the  constant  most  convenient  to  eliminate®  This  can  be  done 
by  dividing  through  by  Ago  Rearranging  and  placing  a  bar  over  the  remaining 
constants  to  indicate  division  by  Ag*  the  equation  becomes 

Bq dT  +  ff  D*~  A-jD©  -  A  ©  -  D2©  (21) 

Substituting  D  *  io>and  writing  in  matrix  notation  gives 


The  normalized  equations  for  m  data  points  can  be  written  directly  by  the 
method  given  previously® 
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C</  2  £iu*f"ef  £-iw©°/  ^~0»eT 

t  </°W£ (i*^)2  o  £~n4‘i^ 

£  -/  «lwe  £(-iuJ&)  £_©‘(-iuii 

,  r\ 


U-<f 


£(-i«*J§)  £_©“  (-iui©)| 

£:~iw©>(4)  £(4) 2  j 


£■=4©  “d" 

£~^©»iu\/ 

£f*?2©"iw© 

£k?2©'>© 

Expanding  the  scalar  products  under  the  summation  symbols  and  letting 
(f  “  A,  +  IB,-  *  1,  the  normalized  equations  become 

r-  ®  ®  l 


m  0  £u?B^  £— A 

0  £l*>2  £  - u2k •  £-*# 

£fa>B£  £.U?2©2  0 

£~  £-*£*  0  £  ^ 


■2”  A^ 
o 

£U;242 


Similarly,,  the  normalized  equations  for  the  constants  in  the  angle  of 
attack  and  normal  acceleration  transfer  functions  are  found  to  be 


m  0  ^  k 

0  2w2  £  -4a  ^“Boc  Cp 

£uJEjjjr  £-Wf^gc  £w>2  oc  0 


ST~4x  ^-^DC 


^©C2 


£-4^ 


<22 

£-U)<x 


m 

0 

£  ~An 

JTuiBn 

0 

£w2 

£“‘4n 

£-h?An 

0 

£~An 

2 

£  n 

0 

^  (uj2a 

n 

£^Bn 

£~^An 

0 

£ui2n2 

£  -4b 

I 

£  H*? 

0 

£id2An 

£-4Bn 

£  . 

£~u?kn 

£-4 Bn 

£w2n2 

0 

£«4a 
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APPLICATION  OF  METHOD  AND  DISCUSSION  OF  RESULTS 


The  least/-squares  method  for  vectors  was  applied  to  actual  flight  test 
longitudinal  response  data  to  obtain  values  for  the  stability  parameters ,  The 
test  aircraft  was  a  twin-engine  light  bomber- type  B-25J  and  is  pictured  in 
Figure  1,  Appendix  I,  A  complete  set  of  frequency  response  characteristics 
for  the  air  speedy  altitude*  and  center  of  gravity  range  of  this  airplane  has 
been  measured  in  flight.,  These  measured  responses  are  intended  to  provide  a 
store  of  flight-test  data  for  the  development  of  data  reduction  methods  and 
to  investigate  the  possibility  of  extrapolating  the  results  of  limited  test 
programso  The  test  methods  and  instrumentation  are  described  in  Reference  5° 

The  response  data*  used  to  demonstrate  the  application  of  the  methods  of 
this  report*  were  obtained  at  an  altitude  of  10*000  feet*  155  mph  indicated  air 
speed*  and  with  the  center  of  gravity  located  at  27  per  cent  of  the  mean  aero¬ 
dynamic  chord„  The  low  air  speed  was  selected  to  provide  a  large  range  of 
values  of  the  reduced  frequency*  K,  where 
UC 

K  =.  2V  (2?) 

CJ  is  the  frequency  of  oscillation  in  radians  per  second,, 

C  is  the  wing  chord,  feet 

V  is  the  free  stream  velocity*  feet/second  ‘ 

Unsteady  flow  effects  are*  theoretically*  a  function  of  the  reduced  frequency . 

The  results  obtained  from  this  example  will  give  some  indication  of  the  possibili¬ 
ty  of  evaluating  these  unsteady  flow  effects  from  flight-test  data. 

The  measured  frequency  response  characteristics  of  the  aircraft  pitching 
velocity  and  normal  acceleration  to  sinusoidal  elevator  oscillations  are  given 
in  Figures  2  and  3  respectively.  Some  of  the  scatter  in  the  data  is  caused  hy 
the  small  variations  in  air  speed,  altitude,  and  weight  that  occurred  during  the 
course  of  the  oscillation  tests.  The  parameters  and  responses  obtained  at  each 
frequency  are  listed  in  Table  1,  Columns  1,  2,  g,  9s  15s  19.--  and  [|2,, 

Table  1  also  presents  a  procedure  for  obtaining  all  of  the  scalar  products 
and  sums  of  products  required  for  the  normalized  least-squares  equations  for 
the  lift  and  moment  derivatives  and  the  pitching  velocity,  normal,  acceleration, 
and  angle- of-attack  transfer  function  constants.  When  only  the  stability 
derivatives  are  required,  this  table  can  be  reduced  from  columns  to  [|0, 

A  convenient  short-cut  method  of  solving  simultaneous  linear  equations 
on  a  desk  calculator  is  outlined  in  Table  2  and  is  self-explanatory.  The 
•theory  and  development  of  this  method  are  explained  in  References  6  and  7» 
Basically,  ii^  consists  of  combining  the  steps  that  are  normally  required  in 
reducing  the  equations,  one  variable  at  a  time,  into  a  complete  operation  for 
each  element.  The  advantage,  when  using  a  desk  calculator,  is  that  the  results 
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of  intermediate  steps  are  retained  in  the  machine,  thereby  eliminating  writing 
down  the  results  of  each  step,, 

The  solution  of  the  normalized  equations  of  motion  for  the  stability 
derivatives  is  outlined  in  Table  3°  The  augmented  matrix,  in  this  case  con¬ 
sists  of  the  matrix  of  coefficients  together  with  the  right-hand  columns  for 
the  lift  equation  and  the  moment  equation  The  matrix  of  the  coefficients 
will  be  the  same  for  the  lift  and  moment  equations  as  shown  in  the  derivation 
(Appendix  I) »  The  auxiliary  matrix  for  the  coefficients  needs  only  to  be  cal¬ 
culated  once  to  obtain  both  the  lift  and  the  momarb  derivatives,,  3he  second- 
right  hand  column  in  the  auxiliary  matrix  is  found  by  the  same  steps  as  the 
first o 


By  arranging  the  variables  in  the  equation  of  motion  in  the  order  of 
descending  magnitude  before  normalizing,  it  becomes  possible  to  check  the 
effect  of  neglecting  the  smaller  derivatives  on  the  values  obtained  for  the 
remaining  derivatives,.  This  can  be  done  by  assuming  that  the  values  of  the 
derivatives  to  be  neglected  are  zero  and  completing  the  remaining  solution  as 
before,,  This  is  shown  in  Table  3  where  C^o  is  assumed  equal  to  zero  and  values 


for  the  remaining  lift  derivatives  are  found,,  Then  Ct°  and  Ct  ,  are  both  dropped 
to  obtain  a  solution  when  only  Ct  is  left0  This  method  will  not  work  if  inter¬ 
mediate  derivatives  are  assumed  ?iro  and  the  remaining  derivatives  to  the  right 
are  not  zero0 


The  values  of  the  derivatives  listed  in  Table  4  agree  closely  to  values 
obtained  by  other  methods  where  comparable  data  are  available „  The  slope  of 
the  lift  curve,  C^,  was  found  to  be  5„111  compared  to  the  value  of  5„25  from 
wind  tunnel  results  (Reference  2>)  „  The  value  of  the  pitching  mom  ait  deriva¬ 
tive,  Cjj^,  equal  to  ~o553s  agrees  well  with  the  wind  tunnel  value  of  -<>525' 
(Reference  g)  „ 

Theoretically,  the  lift  and  moment  resulting  from  elevator  deflection  are  t 
related  as  follows  (Reference  h) s 

c-v'  -ltXV°  '  ~2‘57GV 

where  C  K  the  mean  aerodynamic  chord  «  9 069  feet  for  the  B-25J 

l_k  =  the  effective  tail  length  ■  25„9  feet  for  the  B-25Jo 

It  can  be  seen  that  the  results  in  Table  3  where  Cm  -  -1  „ hlS  and  *  0„556 
are  in  good  agreement  with  this  equation.  The  same^ relationship  should  hold 
for  the  pitching  velocity  derivatives,  however,  the  agreement  is  not  as  good 

here  for  Cmo  »  -0»270  and  CT  „  *  O0II4I0 
rae  ■*->© 

As  pointed  out  in  the  previous  section,  the  effect  of  the  rate  of  change 
of  the  angle  of  attack  is  included  by  an  assumed  linear  relationship  with  the 
rate  of  pitch  derivatives „  The  derivative,  Ct.s  is  so  small  that  it  can  be 
neglected  without  appreciable  effects  on  the  values  for  the  remaining  lift 
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derivatives  o  This  was  verified  by  actually  carrying  out  the  solution  and  ob¬ 
taining  values  of  5»06l,  0O526,  and  O019<4  for  Ct  }  ,  and  C^o  respectively® 

«x  © 


When  the  lift  caused  by  pitching  velocity  is  not  included,,  the  value 
found  for  ,  is  considerably  reduced  (0®2g7)s  although  the  relative  effect 
on  CT  is  much  less  (5®206)®  The  value  for  the  slope  of  the  lift  curve  when 


the  angle  of  attack  is  the  only  variable  included,  is  5»1?60 


V 


One  of  the  reasons  that  the  value  found  for  Ct  is  only  slightly  affected 
by  neglecting  the  smaller  lift  derivatives  is  that  ^he  sums  of  the  products 
in  the  normalized  equations  are  predominately  composed  of  the  responses  in  the 
low  frequency  range®  At  these  frequencies,  the  angle  of  attack  effect  is  large, 
whereas  the  amount  of  lift  from  the  elevator  deflection  and  the  rate  of  pitch 
are  approximately  equal  and  130  degrees  out  of  phase,  thereby  having  little 
net  effect®  At  higher  frequencies,  these  latter  derivatives  become  increasingly 
important  as  their  associated  variables  are  relatively  large  and  their  relative 
phase  angle  changes® 


At  the  high  frequencies,  the  effects  of  other  higher  order  derivatives 
such  as  the  lift  and  momeit  resulting  from  the  angle  of  attack  acceleration, 

Cj^  and  Cj^„  ,  should  become  appreciable®  To  examine  these  effects  the 

normalized  equations  were  obtained  from  the  results  of  a  small  range  of  frequency 
and  solved  as  in  the  example  in  Table  3o  The  sums  of  the  products  from  six  ad¬ 
jacent  frequency  points  such  as  1  through  6,  2  through  7,  etc®,  were  taken  from 
Table  1®  The  values  of  the  derivatives  found  in  this  way  are  plotted  versus  the 
median  frequency  of  each  group  in  Figures  J4.  and  5®  The  derivatives  found  in 
Table  3  are  plotted  as  dashed  lines  over  the  entire  frequency  range  to  provide  a 
comparison  with  the  overall  averages® 


The  angle  of  attack  stability  derivatives  show  consistent  trends  which  can 
be  attributed  to  the  angle  of  attack  acceleration  derivatives®  The  variables 
a  and  ec  are  vectors  along  the  same  line  but  130  degrees  out  of  phase®  The  re¬ 
sponse  caused  by  the  &  stability  derivatives  will  either  increase  or  decrease 
the  effective  value  of  these  stability  derivatives  depending  on  the  relative 
signs®  The  remaining  parameters  are  surprisingly  constant  throughout  the 
frequency  range,  considering  that  the  theory  of  unsteady  flow  indicates  varia¬ 
tions  of  these  derivatives  with  frequency®.  The  scatter  in  the  low  frequency 
range  is  believed  to  be  caused  bv  errors  in  the  data  resulting  from  slight-  tur¬ 
bulence,  although  variation  in  the  speed  of  the  aircraft  during  the  oscillation 
cycle  may  have  contributed  also® 


The  normalized  lift  equations  for  these  limited  frequency  ranges  were  also 
solved  neglecting  first  the  C^®  andCT<>  derivatives  and  then  also  the  deriva¬ 
tive®  The  results  are  plottedein  Figure  6®  The  importance  of  including  these 
latter  derivatives  in  the  solution  at  high  frequencies  is  clearly  shewn  in  this 
graph®  Values  found  for  the  slope  of  the  lift  curve  are  more  than  double  the 
theoretically  expected  size  in  the  high  frequoicy  range  when  these  derivatives 
are  neglected® 


The  effectiveness  of  the  least-squares  process  in  fitting  the  equations  of 
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motion  to  the  original  data  can  be  examined  by  plotting  the  phase  angles  and 
amplitudes  of  the  residual  vectors  for  all  frequencies 0  If  the  equations  of 
motion  include  all  of  the  forces  that  act  on  the  aircraft,  the  residuals  should 
consist  of  small  vectors  with  random  phase  angles  caused  by  the  normal  dis¬ 
tribution  of  errors  in  measuranent-  If  consistent  trends  do  appear  in  the  plots 
of  the  phase  angles  of  the  residuals  versus  frequency ,  it  should  be  possible  to 
find  another  stability  parameter  for  the  equation  of  motion-  This  parameter 
should  be  nearly  in  phase  with  the  residuals-  This  effect  is  demonstrated  in 
Figures  7  a  and  7  b,  where  the  phase  and  amplitudes  of  the  residuals  are  plotted 
for  several  lift  equations-  The  lift  derivatives  are  obtained  from  Table  3» 


The  top  graphs  in  Figures  7  a  and  7  b  are  the  amplitudes  and  phase  angles 
of  the  residuals,  respectively,  whai  only  the  lift  from  angle  of  attack  is  con¬ 
sidered-  !he  amplitudes  are  quite  large  and  the  phase  angles  are  mostly  in  the 
vicinity  of  200  degrees-  Now,  since  exacts  along  the  0  or.  180  degree  vector 
line.,  it  would  appear  that  a  value  could  be  obtained  for  this  parameter-  The 
effect  of  including  ,  in  the  equation  of  lift  is  shown  by  the  residuals  in  the 
middle  graphs ,  Figures^  a  and  7  b-  The  fit  has  been  improved,  as  shown  by  the 
decrease  in  the  sum  of  squares  of  the  amplitudes-  However,  the  phase  angles 
still  show  a  fairly  consistent  trend-  The  inclusion  of  the  G^,,  and  C^.  deriva¬ 
tives  reduces  the  anplitudes  of  the  residuals  still  further  whfle  the  scatter 
in  the  phase  angles  becomes  considerably  more  random- 


It  is  possible  that  further  analysis  will  reval  other  derivatives  that 
will  improve  the  fit  of  the  lift  or  moment  equations  to  the  data-  Care  would 
have  to  be  exercised  to  avoid  introducing  linear  dependencies  in  the  variables 
of  motion  in  order  that  valid  results  would  be  obtained- 

The  effects  of  some  consistent  errors  in  the  measured  flight-test  data 
can  be  pointed  out  at  this  time-  Calibration  errors  in  air  speed,,  altitude, 
temperature,  weight,  and  the  moment  of  inertia  will  be  consistent  for  all  fre¬ 
quency  points-  These  items  appear  primarily  on  the  right-hand  side  of  the 
normalized  equations  of  motion-  ^ince  the  equations  are  linear  these  errors 
will  affect  all  of  the  stability  derivatives  by  the  same  percentage-  The  true 
air  speed  is  also  involved  in  the  computation  of  the  angle  of  attack,  so, 
additional  error  will  result  if  this  item  is  not  accurate-  Control  surface  de¬ 
flections  are  difficult  to  measure  accurately  because  of  twist  and  distortions 
in  the  surface  caused  by  the  air  loads-  Since  only  the  amplitude  is  usually 
in  error,  the  effect  on  the  accuracy  of  the  derivatives  is  minor-  ^he  surface 
effectiveness  will  be  in  error  inversely  as  the  error  in  measurement  while  the 
other  parameters  will  not  be  affected- 


The  transfer  function  constants  were  also  computed  from  flight  test  data 
using  average  values  of  the  flight  conditions  and  flight  test  values  of  the 
stability  derivatives-  These  computations  are  shown  in  Table  7»  ^he  computed 
values  of  the  constants  are  compared  in  Table  8  with  those  obtained  directly 
from  the  normalized  equations-  The  agreemait  is  good  for  most  constants  al¬ 
though  there  is  appreciable  difference  in  several  cases  such  as  Bg-  The  effect 

of  these  differences  was  examined  by  plotting  the  transfer  functions  for  the 
pitching  velocity  in  Figure  8-  The  results  are  in  veiy  close  agreement  throughout 


lh 
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the  frequency  range.  The  actual  flight  test  results  are  plotted  in  Figure  g 
also  to  show  the  *fitw  obtained  by  the  two  methods.  The  transfer  function 
equations  fair  the  .flight  test  data  very  closely  in  all  respects  except  for 
the  amplitude  ratio  in  the  low  frequency  range. 

The  individual  data  points  could  be  weighted,  prior  to  normalizing,  by  a 
factor  proportional  to  their  relative  accuracy  of  measurement.  It  can  be 
seen  that  this  would  improve  the  accuracy  of  the  parameters  that  are  obtained 
in  the  solution  since  the  effect  of  the  less  accurate  data  in  the  normalized 
equations  would  be  decreased.  The  errors  in  the  data  obtained  from  flight 
test  are  predominantly  of  an  incremental  nature  resulting  from  reading  in¬ 
accuracies,  In  the  normalized  equations  for  the  stability  derivatives,  the 
weighting  is  automatically  included  by  summing  the  data  directly  as  measured. 

In  the  normalized  equations  for  the  transfer  function  constants  the  results 
of  high  frequency  oscillations  comprise  the  largest  part  of  the  coefficients  as 
a  result  of  the  high  powers  of'  the  frequency  that  are  involved.  Since  the  re¬ 
sponses  in  this  frequency  range  are  small  these  results  are  probably  less  ac¬ 
curate  and  therefore  should  have  less  weight.  That  this  might  be  the  case  is 
indicated  by  the  comparison  of  the  plots  of  the  transfer  function  equation 
using  constants  from  the  least- squares  solution  and  of  the  actual  .flight  data 
(Figure  g) , 

The  effect  of  increasing  the  weight  of  low  frequency  data  points  was  in¬ 
vestigated  to  determine  whether  the  wfitw  in  this  region  could  be  improved. 

This  was  done  quite  simply  by  including  only  uhe  first  17  data  points  in  the 
normalized  equations.  These  results  are  also  included  in  Table  g.  In  most 
cases  these  values  are  in  close  agreement  with  those  obtained  from  the  previous 
solution.  The  frequency  response  characteristics  for  the  pitching  velocity  from 
these  constants  corresponded  so  closely  to  those  already  shown  in  Figure  g  that 
it  was  not  possible  to  plot  the  results  as  a  separate  curve. 

The  consistent,  results  obtained  throughout  this  investigation  give  a  strong 
indication  of  the  validity  and  accuracy  of  the  aerodynamic  theoiy,  flight  test 
measurements,  and  data  analysis  methods.  For  example,  the  constants,  and  » 
are  common  to  each  of  the  three  transfer  function  equations  investigated.  The 
data  for  the  rate  of  pitch  and  the  normal  acceleration  were  obtained  simultan¬ 
eously  but  by  independent  measurements.  Yet  the  agreement  in  the  values  of  thes 
constants  obtained  from  these  two  sources  is  quite  good,  as  also  is  the  value 
found  by  substituting  stability  derivatives  and  flight  conditions  into  the  theo¬ 
retical  expression  for  these  constants  (Table  S) , 

The  plots  of  the  derivatives  in  Figures  4  and  5  are  other  illustrations  of 
these  consistent  results.  Since  the  data  selected  for  each  set  of  derivatives 
covered  only  a  very  small  frequency  range  the  equations  for  the  six  data  points 
in  each  group  are  nearly  the  same  in  magnitude.  In  other  words,  the  equations 
are  nearly  parallel  and  therefore  the  results  are  highly  sensitive  to  errors. 

The  values  plotted  in  these  figures  show  little  scatter,  however.  The  magni¬ 
tude  of  the  normal  acceleration  residuals  in  the  lower  graph.  Figure  7  a,  com¬ 
pared  to  the  size  of  the  measured  data  in  Figure  3  is  another  example.  Note 
that  the  scale  in  Figure  7  a  is  five  times  as  large  as  Figure  3<> 
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The  method  of  least  squares  for  vectors  also  provides  a  possibility  of 
solving  for  the  aircraft  longitudinal  moment  of  inertia0  The  transfer  function 
Ag  is  directly  proportional  to  the  moment  of  inertia0  All  the  remaining  con¬ 
stants  are  divided  by  Ap  in  the  solution  of  the  normalized  equations »  There¬ 
fore,  the  values  found  for  those  constants  which  are  not  a  function  of  Iy  before 
division  by  Ap  'will  be  inversely  proportional  to  the  moment  of  inertia0  These 
include  X  ,  B^,  I,  C0  20s  and  By  keeping  all  the  other  flight  conditions 
constant  and  obtaining  data  at  two  values  of  the  moment  of  inertia  differing 
by  a  known  incranent,  it  is  possible  to  arrive  at  an  explicit  solution  for  the 
total  moment  of  inertia  for  each  of  these  constants „ 
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CONCLUSIONS 


The  method  of  least-squares  for  vectors  provides  a  method  for  evaluating 
aircraft  longitudinal  stability  derivatives  from  frequency  response  data  pro¬ 
vided  the  linear  dependency  existing  between  the  variables  of  motion  is  eliminated,. 
The  assumption  of  a  linear  relationship  between  the  angle  of  attack  rate  and 
pitching  rate  lift  and  moment  derivatives  eliminated  this. dependency  anc^  gave 
acceptable  results „  This  method  of  least-squares  can  also  be  used  to  find  values 
from  flight  test  frequency  response  data  for  the  constants  in  the  aircraft  tran: 
fer  function  equations 0 

The  lift  caused  by  pitching  velocity  must  be  included  for  conventional  air¬ 
craft  in  order  to  obtain  reasonable  values  for  the  lift  resulting  from  elevator 
deflection „ 

The  use  of  vector  algebra  greatly  simplifies  the  numerical  computation  re¬ 
quired  to  obtain  the  normalized  equations 0 
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APPENDIX  I 


DERIVATION  OF  LONGITUDINAL  DYNAMIC  RESPONSE  EQUATIONS 

The  motion  of  an  aircraft  in  flight  is  prescribed  by  aero  dynamic ,  pro¬ 
pulsive,  gravitational,  and  inertial  forces.  The  equations  of  motion  of  an 
aircraft  for  small  perturbations  from  a  steady-state  flight  condition  have 
been  derived  in  considerable  detail  in  Reference  9«  For  the  dynamic  motion 
considered  in  this  report,  ihe  incremental  changes  in  forward  velocity  are 
considered  negligible,.  The  system  of  axes  are  fixed  in  the  bo<3y  during  the 
motion  with  the  longitudinal  or  X  axis  being  parallel  to  the  flight  path 
during  level  flight.  Nose-up  moments  and  angular  displacements,  velocities, 
and  accelerations  are  positive.  Normal  acceleration  and  force  components  are 
positive  downwards.  Elevator  displacements  and  hinge  moments  are  positive 
when  the  trailing  edge  is  moved  downwards. 


The  equations  for  the  pitching  and  normal  acceleration  degrees  of  free-- 
dom  can  now  be  derived  in  the  form  that  was  found  most  convenient  fop  analyzing 
flight-test  response  data. 


M  -  ije  (l) 

Fj  -  (2) 


where  M  is  the  sum  of  the  aerodynamic  moments  about  the  lateral  axis,  pound 
feet. 

2 

Iy.  is  the  aircraft  moment  of  inertia  about  the  lateral  axis,,  slug-feet 
©  is  the  aircraft  pitching  acceleration,  radian/second2 
Fz  is  the  sum  of  the  aerodynamic  forces  along  the  vertical  axis,  pounds 
m  is  the  aircraft  mass,  slugs 

o 

az  is  the  aircraft  vertical  acceleration,  feet/second 

The  important  incremaitAl  variables  of  aircraft  motion  that  are  capable 
of  producing  aerodynamic  moments  about  the  lateral  axis  are  the  angle  of  attack 
(a),  the  elevator  deflection  (</),  the  pitching  velocity  (§)s  and  the  rate  of 
change  of  angle  of  attack  («)  .  Usually,  only  the  an gl.e-of -attack  and  the 
elevator  deflection  are  considered  important  in  producing  normal  accelerations. 
In  the  higher  frequencies  of  the  range  being  considered,  the  lift  resulting 
from  pitching  velocity  also  becomes  an  appreciable  part  of  the  total  lift  and 
is,  therefore,  included  in  the  development  of  the  equations  of  motion.  The 
lift  resulting  from  a  has  small  effect  on  the  aircraft  motion  but,  since  its 
inclusion  simplifies  the  analysis  of  the  flight-test  data,  it  will  be  retained 
here. 


Equations  1  and  2  can  now  be  written  in  the  form  of  moment  and  lift 
coefficient  equations. 


lg 
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CjJ/’fscec  *  c^JfVW  *  c^JVsce  +  c^/^so*  -  iye 


’  -  -1-2' 

i  =  «:,</,  ©,  or  <sc  are  the  respective  moment  or  lift  co= 


0,  ifv2^  *  CT  JjVW  *  C,  .if^SlS  *  CT  .Ifv2^  -  -mgn 
where  Cm  or  C^>s 

efficients  or,  more  generally,  stabilitv  derivatives 0 
f  is  the  air  density,  slugs/feet^ 

V  is  the  true  velocity,  feet/second 
S  is  the  wing  area,  square  feet 
C  is  the  mean  aerodynamic  chord  (mac) ,  feet 


(3) 
( h ) 


n  is  the  normal  acceleration,  g5s 

g  is  the  acceleration  due  to  gravity,  32°2  feet/second' 


Transposing  the  common  terms  to  the  right  hand  side  and  letting 


h  -  2Iy//V^SC,  second^ 

(5) 

CL  -  2mg//>V2S 

(?) 

the  equations  become 

Sr“  +  °-/+  \  ®  +  Si* 

-  h© 

(7) 

(8) 

If  the  elevator  deflection  <f9  is  considered  as  an  input  forcing  function, 
a  simultaneous  algebraic  solution  of  these  equations  of  motion  will  define  the 
response  of  any  variable  of  the  aircraft  motion,,  This  solution  will  be  ex¬ 
pressed  as  the  ratio  of  two  polynomials  in  D,  the  differential  operator,  d/dt, 
and  corresponds  to  the  transfer  function  equation  as  given  in  servomechanism 
theory o 


The  normal  acceleration,  a  3  is  related  to  the  rate  of  change  in  the 
angle  between  the  horizontal  ana  the  tangent  to  the  flight  path,  §'  „  This 
is  a  result  of  the  fact  that  velocity  is  a  vector  and  an  acceleration  occurs 
when  the  direction  of  the  velocity  is  changed, 

V 

az  *  ng  *  Vif  =  V(dc  =*  0)  or  n  *  g-(dc  -  (9) 

The  equations  can  now  be  expressed  in  operational  form, 

<\  *  V°  * +  (<x ' hD)®  '  -°*/ 

(Glec  +  +  K  +  (CL$  ~ 
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(10) 


where  <7*=  toff  VS,  second,. 

Solution  of  the  operational  equations  by  use  of  determinants  yields  the  transfer 
functions 


© 

B A  +  B,D 

0  1 

" - p— 

(12) 

7 

Ao  +  AiD 

+  k£T 

>* 

(X 

cn  +  C,D 

u  1 

* - - H— ' 

(13) 

c f 

Ao  +  AiD 

+ 

where  AQ 

■  w  W 

+  2T'C 

“k 

m 

A1 

‘  Cte°"‘e 

+  2Ttm  - 

=  hCT  -  C  CT„  +  2fC 

W  %  l©  "k 

(15) 

A2 

-  -hCT  - 
*tk 

■  2Th 

(16) 

B0 

'  V* 

-  CT  c 

(17) 

B1 

=  C  C_ 

<X  cf 

-  CT  G  - 

-  27t 

(lg) 

C0 

-  Cm/L6 

- 

-  276 

V 

(19) 

C1 

(20) 

The  transfer  function  for  the  normal  acceleration  is  obtained  by  combining 


equations 

12  and  1J>0 

s 

V  Ebc  4 

g  /  / 

Eq  +  E^D  +  l^D2 

*  2  ’ 

A0  +  A1B  +  A2D 

(21) 

where  EQ 

*  -VB0/g 

(22) 

h 

-  V(C0  “ 

v*  -  «yL*  -  vv  ■  cv%  *  Cl*V/e 

(23) 

h 

-  TCj/g 

“  VhCT  /g 

V 

(21+) 

The  amplitude  ratio  and  phase  angle  lag  of  any  variable  in  response  to  a 
sinusoidal  oscillation  of  the  elevator  can  be  obtained  from  the  transfer 
function  equations „  These  will  be  derived  for  the  pitching  rate  only,  as  the 
procedure  is  the  same  for  the  remaining  variables  „  The  elevator  angle  is 
expressed  as  a  function  of  time 

/-  I/I  sin  Wt  (25) 

where  fj'j is  the  peak  positive  value  of  the  elevator  deflection  in  radians.. 
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(26) 


LJ  is  the  angular  velocity  of  the  oscillation  in  radians  per  second* 

Then  ©  =  Bn  +  B1D  g  l/|  sin  u>t 

A  +  At  D  +  AoB 
0  J- 


For  sinusoidal  motion  the  differential  operators  D  can  be  replaced  by  i ufs 
where  i  is  the  imaginary  unit,,  ~pro 


Now 


Bq  + 

(Aq  “=  A 2^*^)  +  iWA-^ 


| /|  sinwt 


(27) 


The  complex  numbers  in  the  numerator  and  denominator  can  be  expressed  as 
vectors  in  the  complex  plane  having  an  amplitude  and  a  phase  angle,  that  is. 


A  +  iB  « 


The  condition  that  must  be  satisfied  in  order  that  a  unique  set  of  values 
will  exist  for  the  stability  derivatives  in  equations  7  and  g  is  that  the 
coefficients  of  the  derivatives  must  be  linearly  independent,.  This  is  expressed 
mathematically  by  stating  that  the  determinant  of  the  coefficients  of  n  equa¬ 
tions  with  n  unknowns  shall  be  non-zero 0  The  determinant  of  coefficients  that 
must  be  considered  in  this  case  is  the  same  for  the  moment  and  lift  equations 
and  is  , 

1  ©l  *1 

&~2  d  2  ^2  2 

GCj 

Qb*b 

This  determinant  can  be  shown  to  be  equal  to  zero  because  of  the  relationship 
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given  by  equations  g  and  9°  Solving  equation  g  foroc  gives 


CL 

ni 

V 


°V/i  Ch  z 

ec  "a 


« 


(31) 


where  = 


^lni  +  *^i  +  ^3®i  +  ^If^i 
CL 


etc* 


lj  2$  3j  I4.0 


Substitution  of  equation  31  in  A  gives 


A 


(KiniVi©i«i 

(K3©i)/iei«i 

(KifiUi®A 

( Kin2V2d2«2 

+ 

(  k2<^2)</  2e2^2 

+ 

(K3  ©2)^292*2 

4- 

(Kif  2)<f 2^2*2 

(Kin3^®^3 

(k2/3V3^3*3 

(k3§3V3©3^3 

(Kif3  V3&3 

(K 

(32) 


The  three  determinants  to  the  right  are  equal  to  zero  since  each  has  two 
columns  that  are  in  constant  proportion 0  In  the  same  manner  the  remaining 
determinant  can  also  be  shown  to  be  equal  to  zero  by  the  substitution  of 

1(^-4)  for  n^  in  equation  32° 

g 
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TABLE  1  CALCULATION  AND  SUMMATION  OF  SCALAR  PRODUCTS  FOR  THE  NORMALIZED  EQUATIONS 
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TABLE  3 

SAMPLE  SOLUTION  OF  NORMALIZED  EQUATIONS  FOR  STABILITY  DERIVATIVES 


Columns  Summed  in  Table  1  to  Obtain  Augmented  Matrix; 


0 

(§) 

©  -© 

-© 

0 

m 

©)  - 

-© 

® 

@  -© 

© 

(See  Theory,  Equation  19) 

Augmented  Matrix; 

■1 

20.72335 
-2 „ 13651+ 
•17.89213 

-2o1065l+ 

22 

-1+3.35781+ 

17.89213  107.27231 
-43.35784  -4.65182 
I69. 81986  91.31242 

-13.30576 

-13.33820 

5.79131 

Auxiliary  Matrix; 

20 „ 72335 
-2.10651+ 
17.89213 

-.10165 

21.73587 

-1+1 .53912 

.86337  5.17634 
-I.90670  .23699 
75.16967  .14125 

-.64206 

-.90383 

26959 

Solutions  from  Auxiliary  Matrix 

%  ' 

.1I+125 

°L  * 

.23699  -  .11+125  x 

5.17634  -  0 14125  x 

(-1.90670)  =  .55631 

.36337  -  .55631  x  (-.10165) 

«  5.11094 

cc 

- 

-.26959 

CV  = 

-.90383  -  (-.26959) 

x  (-1.90670)  =  -1.41786 

C  * 

% 

-0 64206  -  (-.26959) 

x  .36337  -  (-1.41736)  x  (- 

.10165)  -  -. 

If  CL. 

is  assumed  to  =  0,  then; 

e 


CT  -  .23699 

C.  -  5.17634  -  .2S699  x  (-.10165)  «  5.20551 

(be 

If  CT  and  C,  are  assumed  to  =  0,  then; 

L©  V 

c  *  5.17634 
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TABLE  k 

SAMPLE  SOLUTION  FOR  THE  PITCHING  VELOCITY  TRANSFER  FUNCTION  CONSTANTS 


Columns  Summed  in  Table  1  to  Obtain  Augmented  Matrix: 


m 

0 

© 

© 

-(2d 

0 

© 

® 

-(3) 

© 

-© 

© 

0 

0 

-© 

0 

© 

© 

(Reference  Theory,  Equation  24) 

Augmented  Matrix; 

22 

0 

100.54611 

30.15334 

295.41151 

0 

679.0712 

295.41151 

- 

100.54611 

-4690.1187 

100.51+611 

295.41151 

960.15451 

0 

0 

30.153^ 

-100.514.611 

0 

88.81077 

960.15451 

Auxiliary  Matrix: 

22 

0 

4.57027 

1.37062 

13.42779 

0 

679.0712 

.43502 

-.11+806 

-6.90666 

100.51+611 

295.41151 

372.12172 

-.25279 

1.35475 

30.1538I+ 

-lx.  54611 

-94.07162 

8.811+09 

4.00456 

Solutions  from  Auxiliary  Matrix: 

AQ  =  1*.  00456 

I  =  1.85475  -  4.00456  x  (-,25279)  =  2.S6706 

B  -  -6.90666  -  4.00456  x  (-.114806)  -  2.86706  x  .1+3502  «  -7.56097 

Bq  -  13.142779  -  J+.Xl;56  X  1.37062  -  2.86706  X  U*  57027  -  -5.161+13 
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TABLE  5 


SAMPLE  SOLUTION  FOR  ANGLE  OF  ATTACK  TRANSFER  FUNCHON  CONSTANTS 


Columns  Summed  in  Table  1  to  Obtain  Augmented  Matrix: 


m 

0 

-© 

-Q 

-© 

0 

© 

© 

-©  , 

-© 

© 

0 

0 

-dj) 

© 

0 

© 

(See  Theory ?  Equation  25) 

Augmented  Matrix: 

22 

0 

29.314221+ 

2ol065l+ 

-72.87690 

0 

679.0712 

-72.87690 

-29.3I422I+ 

1411.5141432 

29-314221+ 

-72.87690 

73.66311 

0 

0 

2.10651+ 

-29,314221+ 

0 

20.72335 

73.66311 

Auxiliary  Matrix; 

22 

0 

1.33373 

.09575 

-3.31253 

0 

679.0712 

-.10731 

-0  01+320 

-,60601+ 

29  0^14221+ 

-72.87690 

26.70806 

-.22307 

1.93562 

2.10651+ 

-29.314221+ 

-5.95327 

17.925145 

1+.  16667 

Solutions  from  Auxiliary  Matrix; 

*0  * 

1+.  16667 

h  " 

1,93562  -  ( 

-.22307)  x  i+.  16667  =  2,91508 

h  " 

-.60601+  -  ( 

-.01+320)  x  1+.  16667  -  (-.10731) 

x  2.91508  *  -.11322 

C n  " 

-3.31253  -  . 

09575  X  i+.  16667 

-  1.33373  X  2.9150s  =  -7.599I47 
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TABLE  6 


SAMPLE  SOLUTION  FDR  NORMAL  ACCELERATE ON  TRANSFER  FUNCTION  CONSTANTS 


Columns  Summed  in  Table  1  to  Obtain  Augmented  Matrix; 


m 

0 

-© 

© 

-© 

-© 

0 

© 

-© 

0 

-© 

-© 

-© 

© 

0 

© 

® 

© 

~© 

0 

© 

-(§) 

0 

-© 

0 

© 

© 

® 

(Reference 

Theory,  Equation  26) 

Augmented  Matrix; 

22 

0 

-6.91J+52 

-227,21209 

-679^0712 

0 

679,0712 
227, 21S 09 
955,1+7050 
0 

-6,911452 

227,21209 

1070,171*51+ 

0 

-955,1+7050 

-227,21209 

955,1+7050 

0 

1+91+1,9625 

1+195,3905 

-679,0712 

0 

-955,1+7050 

1*195,3905 

1*2370,050 

955,1+7050 

1*195,2905 

1*91+1,9625 

0 

-63366,252 

Auxiliary  Matrix; 

22 

0 

-6,911+52 

-227,21209 

-697o0712 

0 

679,0712 

227o21S09 

955,1+7050 

0 

-,311+29 

,331+60 

991,971+20 

-391,11220 

-1162,29579 

-10,52209 

1,1*0702 

-,391+27 

109606631*7 

-3273,1+7351 

-30,26627 

0 

-1,17335 

-2,92951 

16230,231*9 

1+3,1+301+7 
6,17226 
3,26936 
1+,  991+99 
-,7991+0 

Solutions  from  Auxiliary  Matrix; 


E2  =  7991+0 

A  =  2 o6051g 

A  =  3»95l+53 

0 

E1  -  lo 19013 

Eq  »  1+6  ,,9014.90 
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TABLE  7 


TRANSFER  FUNCTION  CONSTANTS  COMPUTED 
FROM  FLIGHT  DATA 


Parameters  Obtained  From  Flight  Test  Datas 


h 

ave 

2 

*  -ISO  sec 

11 

6« 

0 

”-555 

V 

5-111 

T  ave 

*  3-00  s  ec 

c 

— I0I4IS 

V 

-556 

V 

ave 

»  265  ft/sec 

c 

m° 

© 

=>-270  sec 

V 

0II4O  sec 

g 

-  32 o2  ft/sec2 

C 

m  - 

oc 

~-122  sec 

cfe ' 

-063  sec 

Transfer  Function  Constants  Calculated  from  Flight  Test  Parameters* 


So 

BX 

*0 

^1 

*L 


-hCT 

L° 

OC 


-1-091  secJ 


<C1  Cm.  -  \  *  2Tcm  )/k 

OC  ©  QC  •  <*  c 


(CT  C,,,  +2  TCm 


(CT0 
L ,  ® 
<f  oc 


sec 

hC t  -  C  CT  +  2Tc  /A0  »  2<>99g  sec 
L  mo  L”  2 

oc  oc  ®  oc 

CT  C  )/A„  «  -6-359  sec”5 


(c„cT  -cTac. 


)/A  *  “7oSl6  sec 


-2 


»  -7-752  sec" 


V  ¥  m<f  ’■  m<f  2 

(C  C.  *  C.C  -  2TCm  )/Ap 

hC,  /A  «  =0-092  gee'1 

V  2 

V(SCW  *  52.33U  =eo-2 

7 ^ C  CT  -  Ct,  C„0  +  Ct  C_  +  Ct  „C_  ) /eA„ 

VhCT  /gA0  «  -Oo755 


0-527  sec 
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TABLE  2 


TRANSFER  FUNC  HON  CONSTANTS  FROM  FLIGHT  TEST  DATA 


(Bq  +  B1D)/(I0  + 


XxD  +  D2) 


B, 


B, 


_0 


From  Data  Points 
1  through  22 

bo005 

2.267 

-5.16U 

-7.561 

(co  ♦  ciD)/(I0  +  V  ♦ 

Ud67 

2.915 

“7.599 

-0.113 

d0  +  E-jD  +  &/)/(AQ 

3«95i+ 

2.605 

U6.905 

1.190 

“0.799 


From  Data  Points 
1  through  17 

3.3140 

2.917 

-5.116 

-7.772 

>2) 

I4.I7U 

2.959 

-7.732 

-0.122 

2. 

A-jD  +  D  ) 

U.122 

2.229 

51.211 

0.1^25 

-0.655 


From  Stability- 
Derivatives  ,  Table  7 

I4.23I4 

2.99^ 

=6.399 

-7.216 


“7.752 

-0.092 


52.3314 

0.527 

-0.755 
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